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We illustrate how attosecond light pulses can be used to directly mapping out the time-dependence 
of the correlated motion of two excited atomic electrons, discuss how the two-electron correlations 
manifest themselves in realistic attosecond measurements, and propose the following for experimen- 
tal exploration: (a) The single ionization signals which directly reveal bending-vibrational motion 
of the correlated electron pair, (b) and also its rotational motion, (c) The double ionization signals 
which directly reveal the two-electron density in momentum space. To facilitate the description of 
the above points, use is made of simple wave packets of doubly-excited states of helium. 
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Following the recent developments in intense laser 
physics, ultra-short XUV (extreme ultraviolet) light 
pulses with duration of several hundred attoseconds (or 
a few atomic units) have been reported 0, 0. Such 
pulse durations are comparable to the time scale of the 
electronic motion in the ground and the lower excited 
states of atoms and molecules, thus opening up the route 
to the time-resolved study of electron dynamics in mat- 
ter, akin to the time-resolved tracking of the atomic mo- 
tion in a molecule enabled by the advent of femtosec- 
ond laser pulses ^4]. In terms of probing electron dynam- 
ics, the only time-domain measurement reported so far is 
the determination of Auger lifetime by Drescher et al. , 
although other applications of attosecond pulses have 
been reported elsewhere |^ 0. Unlike molecules where 
the time-dependent rotational and vibrational motions 
have clear classical meaning, the significance of the time- 
dependence of the electronic motion is presently elusive. 
The root of this misapprehension lies in the shell model of 
atoms which currently forms the basis for interpretation 
of almost all the energy-domain measurements. Accord- 
ing to this model, each electron is moving in an effective 
central-field potential made of the electron-nucleus in- 
teraction plus an average potential due to the remaining 
electrons. Thus the time dependence of the interaction of 
one electron with the others is often probed in the form 
of a relaxation, which tends to be monotonic. Standard 
atomic structure theory accounts for the deviations of the 
central field model in terms of configuration interaction 
(CI). Thus in the time-domain study of electron dynam- 
ics, the time-dependent CI coefficients are calculated ^T^. 
The point is that such coefficients carry little physical 
meaning, and one must resort to some alternative view- 
point for a physical interpretation. 

In this Letter we suggest that a good place to study 
the time-domain electron dynamics is where the break- 
down of the shell model is most severe, namely to probe 
multiply excited states of an atom using attosecond XUV 
pulses. Theoretical studies in the past decades have re- 
vealed that the shell model fails completely for these 



states H i, [T3, 111 111 El El El El, and the motions 
of electrons in these states are better described by draw- 
ing analogy with the rotation and vibration of a floppy 
polyatomic molecule, where the periods are of the or- 
der of a few to sub- femtoseconds. By using doubly ex- 
cited states of helium as an example, we show that the 
rotational and bending vibrational motions of the two 
electrons can be probed directly with XUV attosecond 
pulses, where single and/or double ionization signals of 
coherently populated doubly excited states should reveal 
these rotational and/or vibrational motions directly. For 
these systems, the "movies" of the correlated motion of 
the two electrons can be made to reveal their time evo- 
lution, similar to the "movies" showing the motion of 
atoms in a molecule. However, there is a signiflcant dif- 
ference. The rotational and the vibrational periods in 
a molecule are at least two orders of magnitude differ- 
ent. Thus their motions are not directly coupled. For 
atoms, the rotational and vibrational periods are compa- 
rable. To disentangle the various rotational and bending 
vibrational modes would require a careful preparation of 
the initial coherent double excited states. Addressing the 
preparation of such coherent doubly excited states is not 
a purpose of this Letter. Instead, we consider simplest 
initial coherent states and calculate the time-dependence 
of the ionization yield to show that they indeed reveal the 
rotational and the bending motion of the two electrons 
directly. 

To describe the collective motion of the two electrons 
in the doubly excited states EH, it is convenient to look at 
them as a linear XY2 triatomic molecule, with X play- 
ing the role of the nucleus and Y an electron. Instead 
of the independent electron coordinates ri and r2, the 
wavefunctions are to be expressed in terms of = 
(a, ^12) for the internal stretching and bending vibra- 
tional motions, and three Euler angles Qr = (<^^/^^70 
for the overall rotational motions. Here the hyperradius 
R and the hyperangle a are deflned hy ri = R cos a and 
r2 = i^sina. Thus R stands for the size of the atom, a 
measures the relative distances of the two electrons from 
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the nucleus, and 612 is the angle which the two electrons 
make with the nucleus at the vertex. Note that each 
wavefunction depends on six variables excluding spins. 
For visualization, the vibrational and rotational mode of 
each state can be defined via the densities. 
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for each state (pj. In presenting we choose ri — 

1*2 and ri x 1*2 to be parallel to the z'- and ?/'-axes of 
the body- fixed frame, respectively. We consider linearly 
polarized laser pulses parallel to the 2:-axis of the space- 
fixed frame. Thus, we only need to consider the magnetic 
component of the total angular momentum, M = 0, so 
that the density does not depend on 7'. 
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FIG. 1: (color online). Density plots of the doubly excited 
states. The top row presents the relief plots of p^^^(a,^i2), 
and the middle row shows p^^^ {a , f3') by polar plots. As seen 
from the density plots, 2s^ ^5'^ 2s^ ^D"", and 2/ ^5'^ are 
considered as ground, rotationally excited, and vibrationally 
excited states, respectively. The bottom row sketches the cor- 
responding classical rotational and bending vibrational modes 
of the two correlated electrons. 

Neglecting autoionization for the time being, consider 
first the "stationary" doubly excited states, convention- 
ally called 2^2 ^S'^ 2^^ ^S"", and 2p'^ states. These 
states exhibit large configuration mixing and are more ac- 
curately classified using the (i^, T)^ quantum numbers 
as described in 11]. In this designation, the three states 
are denoted as {K,T)^ = (1,0)+, (-1,0) + , and (1,0)+, 
respectively. In Fig.^the vibrational and rotational den- 
sities of each state are given, together with a classical 
sketch of the motion of the two electrons. 

For 25^ ^iS^, p^^^(a,6>i2) has a maximum at 6^12 = tt 
and p^^^{^}y) is isotropic. Thus, this state is considered 
as the ground state of the rovibrational motion. For 
2j9^ ^D^^ p^^^{a^0i2) has maximum at O12 = tt simi- 
lar to 2s^ -^5^, but p^^^{ftr) has nodes in similar to 



|^2o(/^^ <^0P total angular momentum L = 2 with 

T = 0. That is, this state is rotationally excited. For 
2p^ ^S^^ on the contrary, is isotropic, but there is 
a nodal line in p^^^(a,6>i2) at about 6^12 = 7r/2, so that 
this state is an excited state in the bending-vibrational 
mode. With this understanding, we proceed to discuss 
how to observe the time-dependence of these rotational 
and vibrational motions by first creating coherent states 
by a pump pulse. We will then consider the use of at- 
tosecond light pulses for probing the time dependence of 
the coherent state via the ionization yield. To this end, 
we calculated the ionization probability from such coher- 
ent doubly excited states by applying the time-dependent 
Hyperspherical methodflTj to the single ionization, and 
the first order perturbation theory to the double ioniza- 
tion separately. 

We highlight the following three points in turn. (1) 
Attosecond pulses probing bending vibrational motion. 
When the 2s^ ^5*^ and 2p^ ^S'^ states are coherently 
excited at t = 0, a vibrational wave packet is created. 
We used a Gaussian attosecond pulse of mean energy of 
21.8 eV, duration of 142 asec, and intensity of 3.5 x lO"*^^ 
W/cm^ to ionize this coherent state by stimulated emis- 
sion. The resulting ionization yield of He+(ls)+e in the 
energy region of —2 to —1 a.u., as shown in the top frame 
of Fig. [3 shows the expected oscillations as a function of 
delay time. The oscillation has a period of about 970 
asec, corresponding to the inverse of the energy separa- 
tion between the two states, 27r/A£^. In fact, this oscil- 
lation can be traced directly to the bending vibrational 
motion, as seen clearly from the calculated time evolu- 
tion of the average angle between the two electrons with 
the nucleus at the vertex, {O12). At the bottom of this 
figure, the calculated p^^^(a,^i2) is shown over a period. 
The ionization probability peaks when the two electrons 
are mostly on opposite sides of the nucleus. We note that 
the ionization yield decreases with increasing delay time, 
due to the autoionization of the 2s'^ ^S^ state which has 
a lifetime of about 5.4 fsec. Thus the ionization yield can 
be used to measure the lifetime of the shorter-lived state 
directly without the help of an IR laser ^J. Note that at- 
tosecond pulses are needed in order to see the identifiably 
rapid oscillations in the ionization probabilities. 

(2) Attosecond pulses probing the rotational motion. In 
Fig. [31 we show the ionization probability of IIe+(ls)+e 
in the energy region of —2 to —1 a.u. vs delay time from 
a coherent state made of 2s'^ ^S^ and 2p^ ^D^. In this 
case each state is the ground state of the vibrational mo- 
tion. The time-dependent oscillation is to be traced to 
the rotational motion, since the 2p^ ^D^ state has total 
angular momentum L = 2 with T = 0, i.e., it is rotation- 
ally excited in /3\ The time dependence of the ionization 
probability is shown to follow the oscillation of the Euler 
angle (The average of the deviation from 7r/2.) This 
angle is defined to be the angle between the "molecular 
axis" (the axis between the two electrons) with respect 
to the laser polarization direction. In this case, more 
ionization occurs when the line joining the two electrons 
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Hons in momentum space. The oscillation of the single 
ionization probabilities presented above is generic to any 
two-state system. The oscillation can be interpreted in 
terms of bending vibrational or rotational motions of the 
two electrons only when they are accompanied by such 
careful theoretical analysis as presented here. Without a 
priori information about the nature of the coherent state, 
a blind display of the six-dimensional two-electron wave- 
function would not be able to reveal the nature of the 
two-electron dynamics. For an arbitrarily created coher- 
ent state, an ambiguous analysis of this sort would be 
rather difficult. We thus ask whether it is possible to 
map out the time-dependence of the two-electron motion 
directly. One method we propose is to measure double 
ionization of helium by an attosecond light pulse. 
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FIG. 2: (color online). Delay time dependence of (a) stim- 
ulated single ionization probability from a coherent doubly 
excited state of [ip{2s'^ ^S'^) + 2(^(2/ ^5'^)] /VS of He, (b) 
the average angle between the two electrons with respect to 
the nucleus, and (c)the two electron density distributions over 
one period. 



are parallel to the laser polarization. A more complete 
view of the rotational motion of the wave packet can be 
seen from the bottom frame where the rotational density 
evolution within one period is displayed. Note that the 
oscillation period (2.0 fsec) is inversely proportional to 
the energy separation between the two states. 
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FIG. 3: (color online). Delay time dependence of (a) stim- 
ulated ionization probability from a coherent doubly excited 
state of [ip{2s^ ^5'^) + ip{2p^ ^D^)] /a/2 of He, (b) the aver- 
age of the angle between the two-electron axis with respect 
to the laser polarization (see text), and (c) the rotational dis- 
tributions of the atom in Euler angles. 
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FIG. 4: (color online). Time-dependence of double ioniza- 
tion from a coherent state of [v?(2s^ ^S"") + ip{2p'^ ^5'^)] /V2. 
(a) The angular distribution of the electron density of the 
initial coherent state in momentum space. Shown is the mo- 
mentum density of the second electron with respect to the 
direction of the fixed first electron; (b) The distribution of 
the calculated ionization signal of the second electron with 
respect to the direction of the fixed first electron; (c) The 
two-dimensional two-electron momentum space density dis- 
tributions of the initial coherent state; (d) Double ionization 
signals showing the relative momentum distributions between 
the two ionized electrons, (see text). 



(3) Attosecond pulses probing two-electron wavefunc- 



We calculated the probability of double ionization by 
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an attosecond light pulse with mean energy of 27.2 eV 
and pulse duration of 83 asec. We use a weak light pulse 
so that double ionization probability is calculated by the 
first order perturbation theory. In the actual calculation, 
the velocity gauge was used and the electric field of the 
attosecond pulse was assumed Gaussian. The two free 
electrons in the final states are approximated by prop- 
erly symmetrized products of plane waves. Under this 
approximation, one can formally show that in the limit 
of a very narrow pulse, the ionization probability is pro- 
portional to the two-electron charge density in the mo- 
mentum space of the initial coherent state. 

We have carried out a numerical calculation for the 
simplest coherent state made of 2s^ ^S^ and 2p^ ^S^. The 
double ionization probability P(pi, p2, t) is a function of 
the six-dimensional momenta of the two electrons. Here 
t is the time that elapsed after the pump pulse. This par- 
ticular example of a coherent state represents a bending 
vibrational wave packet. In Fig. ^a) we plot the ion- 
ization rate 5'(pi,p2, t) ex P(pi,p2,t)/ [e • (pi + p2)]^ as 
a function of Op2^ with fixed ejected electron energies of 
two electrons ei = 62 = 2.2 eV and ^pi = 0, over half a 
period at t = 0, T/8, T/4, 3T/8, and T/2 asec (the total 
double ionization probability is expected to oscillate with 
a period of T = 970 asec for this coherent state) Here, 
the angles are measured relative to the laser polarization 
e. In the upper frame, the density plots constructed from 
the momentum space wavefunction are shown. It is clear 
that the angular dependence of the ionization signal re- 
sembles the angular dependence of the momentum space 
electron density; see Fig. ^b). Note that the amplitude 
of the oscillations in the ionization signal is somewhat 
weaker. This is due to the averaging effect from the fi- 
nite pulse duration of 83 asec. (The mean energy of 27.2 
eV was used, since such short attosecond pulses have to 
be generated from the plateau region of the high-order 
harmonics.) 

We define the momentum space hyperspherical coordi- 
nates for the two electrons which enable us to similarly 



define the momentum space counterpart of the vibra- 
tional and rotational densities. In Fig. ^c) and (d) we 
compare the ionization measurement with the momen- 
tum space densities as functions of = arctan(p2/pi) 
and 6pi2 = arccos(pi • P2)j30- (Note that these an- 
gles are not conjugate to the hyperspherical coordinates 
variables in the configuration space, a and O12.) Indeed, 
the ionization signal clearly reproduces the bending vi- 
brational density of the two electrons in the momen- 
tum space. Representing the final state wavefunctions 
by plane waves as in the present example is undoubtedly 
an over-simplification. Improved calculations might show 
slight distortions to the double ionization signals, a task 
left for future exploration. Nevertheless, it is undeni- 
able that the present calculation points out that attosec- 
ond pulses are capable of mapping out the two-electron 
dynamics. Such information cannot be directly derived 
from energy-domain measurements. 

In summary, we have shown that attosecond light 
pulses can be used to probe the correlated motion of 
two excited electrons. For properly created coherent 
states the rotational and/or bending vibrational modes 
of the two excited electrons can be directly mapped by 
the single or double ionization signals vs the time de- 
lay. Extension of the present analysis to multiply excited 
states and many-body systems would similarly 

offer a new perspective that distinguishes itself from the 
established energy-domain measurements. 
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